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ABSTRACT 

The multiple system V505 Sagittarii is composed of at least three stars: a compact eclipsing 
pair and a distant component, which orbit is measured directly using speckle interferometry. In 
order to explain the observed orbit of the third body in V505 Sagittarii and also other observable 
quantities, namely the minima timings of the eclipsing binary and two different radial velocities 
in the spectrum, we thoroughly test a fourth-body hypothesis — a perturbation by a dim, yet- 
unobserved object. We use an N-body numerical integrator to simulate future and past orbital 
evolution of 3 or 4 components in this system. We construct a suitable x 2 metric from all available 
speckle-interferometry, minima-timings and radial- velocity data and we scan a part of a parameter 
space to get at least some of allowed solutions. In principle, we are able to explain all observable 
quantities by a presence of a fourth body, but the resulting likelihood of this hypothesis is very 
low. We also discuss other theoretical explanations of the minima timings variations. Further 
observations of the minima timings during the next decade or high-resolution spectroscopic data 
can significantly constrain the model. 

Subject headings: stars: binaries (including multiple): close - stars: individual: V505 Sgr 

1. Introduction visual magnitude in maximum 6™ 5 (Chambliss et 

al. 1993). In 1985, the V505 Sgr was also rc- 
The star V505 Sagittarii (HD 187949, HR 7571, solved using spec ki e interferometry (McAlister et 
HIP 97849, WDS 19531-1436) is known as an al 19g7a) and several mea surements of this 3rd 
eclipsing binary with a variable period. Spectral component were published since that time. Mayer 
types of its primary and secondary components (1997) attem pted to join the measured times of 
are A2V and G5IV, orbital period is 1.183 and minima with visual orbit and determined a dis- 
tance of the system 102 pc. 

1 On leave from Astronomical Institute, Slovak Academy . . 

of Sciences, 05960 Tatranska Lomnica, Slovakia The 3rd-body orbit With the period of about 
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40 years seemed well justified until about the year 
2000. An abrupt change in more recent data how- 
ever excludes this simple model — it is impossi- 
ble to fit both light-time effect data and the in- 
terferometric trajectory assuming three bodies on 
stable orbits. We thus test a 4th-body hypoth- 
esis: a perturbation by low- mass star (i.e., the 
4th body), which is not yet visible in the speckle- 
interferometry data. Such a fourth body was sus- 
pected already by Chochol et al. (2006) due to 
conspicuous deviations of minimum times from ex- 
pectations. While we consider the 4th-body model 
as the main working hypothesis in this paper, we 
also discuss other possible effects that can produce 
minima timings variations. 

The dataset we have for V505 Sgr is described 
in Section [2] We introduce our dynamical model, 
numerical method, free/dependent parameters 
and \ 2 metric in Section [3] The results of our 
simulations and conclusions are presented in Sec- 
tions H] and [5j 

2. Observational data 

2.1. Speckle interferometry 

The available speckle-interferometry data are 
summarised in Table [TJ Most of them were ex- 
tracted from the Fourth Catalog of Interferomet- 
ric Measurements of Binary Stars (Hartkopf et al. 
2009), but we also added two speckle measure- 
ments from SAO BTA 6 m telescope by E. Mal- 
ogolovets (using a speckle camera and a method 
described in Balega et al. (2002) and Maksimov 
et al. (2009)) and one direct-imaging measure- 
ment, performed at CFHT by S. Rucinski (using 
a method described in Rucinski et al. (2007)). 

We estimated weight factors w and correspond- 
ing uncertainties as er s k y = 0.005 arcsec/w. The 
values of er s k y vary because of different telescopes 
and techniques were used. Note these measure- 
ment errors sometimes cause that some measured 
points are seemingly 'exchanged' — the position 
angle does not revolve monotonically. 

We are aware of a possible 180°-ambiguity in 
the speckle measurements, but V505 Sgr is a lucky 
case: we have one direct measurement by Hippar- 
cos prior to 2000 perihelion passage, and another 
direct- imaging datum after 2000. We thus can be 
sure about the shape of the orbit. 



2.2. Minima timings 

We list recent O — C data for the (1+2) binary 
in Table [3J Only measurements not presented in 
Chambliss et al. (1993) are included in the ta- 
ble, but we use all of them of course. An uncer- 
tainty of a minimum determination is estimated 
to (Tiito = 1 min in most cases, only photographic 
minima and data from Hipparcos were considered 
worse. Epoch and O — C were calculated using the 
cphcmeris: 

Pri.Min. = 2433490.483 + l d 1828688 x E (1) 

Note there is a freedom in period and base min- 
imum determination. When we compare these 
O — C measurements to our simulations we use 
an optimal ephemeris, different from (TTJ). 

2.3. Radial velocities 

We use radial- velocity data from Tomkin (1992), 
Tab. 4, who measured sharp spectral lines in 
the 5580 - 5610 A region and attributed them 
to the 3rd component. The values of w ra d3 
range from —13 to — 9km/s. The width of the 
lines corresponds to rotational velocity about 
«rot3 = (20 ± 5) km/s. 

The uncertainties of the radial-velocity data 
<t iv = 2kms~ 1 were estimated from a scatter of 
the RV measurements close in time and verified 
by a practical test: we computed a synthetic spec- 
trum with the same resolution as Tomkin (1992) 
and fitted the lines in question by a Gaussian func- 
tion. We also checked for possible blends with 
nearby faint lines. 

Wide lines in the V505 Sgr spectrum are at- 
tributed to the components of the eclipsing pair 
(1+2). The binary is tight and in all likelihood 
rotates synchronously, thus the corresponding ro- 
tational Doppler broadening is large (v ro ti+2 = 
(100 ± 10) km/s). The systemic radial velocity of 
the (l+2)-body is v ladl+2 = (1.9 ± 1.4) km/s. 

3. Numerical integrator and x 2 metric 

In order to model orbital evolution of the 
multiple-star system V505 Sgr, namely mutual 
gravitational interactions of all bodies, we use a 
Bulirsch-Stoer A-body numerical integrator from 
the SWIFT package (Levison & Duncan 1994). 
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Table 1: Speckle interferometry data for V505 Sgr, mainly from the 4th Interfcrometric Catalogue (Hartkopf 
et al. 2009). P.A., d denote position angle and angular distance between the central pair (1+2) and the 3rd 
component. Estimated weight factors w and uncertainties cr s ky = 0.005 arcsec/w correspond to the sizes of 
telescopes and techniques, which were used to acquire these measurements (1991.25 and 2005 measurements 
result from direct imaging). 



year 


P.A. / dc 


g d J mas 


weight source 


1985.5150 


189. 


.6 


302 


1 


3.6 m 


1985.8425 


189. 


.8 


311 


1 


3.8 m 


1989.3069 


181. 


.0 


261 


1 


4.0 m 


1990.3445 


176. 


.9 


246 


1 


4.0 m 


1991.2500 


170 




234 


0.6 


Hipparcos 


1991.3903 


173. 


.4 


234 


1 


4.0 m 


1991.5575 


174 




240 


0.4 


2.1 m 


1991.5602 


174 




260 


0.4 


2.1 m 


1991.7124 


173. 


3 


226 


1 


4.0 m 


1992.4497 


171. 


.7 


214 


1 


4.0 m 


1992.6961 


164 




190 


0.4 


2.1 m 


1994.7079 


159. 


.9 


192 


1 


3.8 m 


1995.4398 


152. 


5 


169 


0.6 


2.5 m 


1995.7675 


154. 


.2 


177 


0.3 


2.5 m 


1996.5320 


145. 


.8 


149 


0.3 


2.5 m 


2003.6365 


236. 


.3 


152 


1 


3.5 m 


2005.7948 


218 




183 


0.6 


direct CFHT 


2006.1947 


215. 


.8 


182 


1 


4.0 m 


2007.3306 


212. 


.4 


210 


1 


3.5 m 


2007.4927 


212. 


.0 


212 


1 


6.0 m 


2008.4901 


207. 


.8 


231 


1 


6.0 m 


2009.2662 


204. 


.2 


247.5 


1 


4.0 m 



3 



Table 2: Minima timings for the eclipsing binary (1+2) in V505 Sgr. Epoch and O — C were calculated using 
the ephemeris Pri.Min. = 2433490.483+ 1^1828688 x E. ante denotes assumed standard uncertainty of the 
minimum determination. Only newer minima after Chambliss et al. (1993) are listed. The references are: 
Rovithis-Livaniou & Rovithis (1992), Miiyesseroglu et al. (1996), Ibanoglu et al. (2000), Cook ct al. (2005), 
Chochol et al. (2006), Zasche ct al. (2009). The last 4 measurements are new. 



JD - 2400000 Epoch 0-C/d ante I d source 



48432. 


.4871 


12632. 


.0 


+0 


.0054 





,0007 


R.-L. 


48501. 


,0981 


12690 


.0 


+0 


.0100 





,0021 


Chochol 


48858. 


.3253 


12992. 


.0 


+0 


.0109 


0. 


,0007 


Miiyesseroglu 


51000. 


.4948 


14803. 


.0 


+0 


.0049 


0. 


.0007 


Ibanoglu 


51051. 


.3578 


14846. 


,0 


+0 


.0046 


0. 


.0007 


•>•> 


51057. 


.2724 


14851. 


.0 


+0 


.0048 


0. 


.0007 


» 


51064. 


.3692 


14857. 


.0 


+0 


.0044 





.0007 


ii 


52754. 


.6756 


16286. 


.0 


-0 


.0087 


0. 


.0007 


Chochol 


52843. 


,3891 


16361 


.0 


-0 


.0103 


0. 


.0007 


ii 


53263 


.3029 


16716 


.0 


-0 


.0150 





.0007 


ii 


53525 


.8969 


16938. 


.0 


-0 


.0178 





.0007 


Cook 


53626. 


.4399 


17023. 


.0 


-0 


.0187 


0. 


.0007 


Chochol 


54267. 


.5469 


17565. 


.0 


-0 


.0266 


0. 


.0007 


Zasche 


54267. 


.5472 


17565. 


.0 


-0 


.0263 


0. 


.0007 


ii 


54648. 


.4260 


17887. 


.0 


-0 


.0313 


0. 


.0005 


ii 


54655. 


.5233 


17893. 


.0 


-0 


.0312 


0. 


.0003 


ii 


54658. 


.4817 


17895. 


.5 


-0 


.0299 


0, 


.0005 


ii 


54706. 


.3869 


17936. 


,0 


-0 


.0309 





.0002 


ii 


55027. 


.5302 


18207. 


,5 


-0 


.0365 


0. 


,0018 


Uhlaf 


55049. 


.4152 


18226. 


.0 


-0 


.0346 


0. 


,0002 


ii 


55062. 


.4266 


18237 


.0 


-0 


.0347 





,0011 


Smelcer 


55068. 


.3400 


18242. 


,0 


-0 


.0357 


0. 


,0011 


Uhlaf 
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Our method is quite general — we can model 
classical Keplerian orbits, of course, but also non- 
Keplerian ones (involving 3-body interactions). 
We are able to search for both bound (elliptical) 
and unbound (hyperbolic) trajectories. Free pa- 
rameters of our model are listed in Table [3l Here- 
inafter, we strictly denote individual bodies by 
numbers: 1, 2 (Algol- type pair), 3 (resolved third 
component) and 4 to avoid any confusion. 

Fixed (assumed) parameters are listed in Ta- 
ble ID Masses of the first three components are 
well constrained by photometry and spectroscopy: 
mi = (2.20 ± 0.09) Mq, m 2 = (1.15 ± 0.05) M Q , 
m 3 = (1.2±0.1)M Q (Chambliss et al. 1993, 
Tomkin 1992). We take 777,3 as a free parame- 
ter, thought, because of larger relative uncertainty. 
When we test 3-body configurations, we have sim- 
ply m 4 = 0. 

First, it is often useful to adopt a simplifica- 
tion: 1st and 2nd body can be regarded as a single 
(1+2) body in our dynamical model. The central 
pair (1+2) is so compact (a = 0.033 AU) and the 
distance of other components so large, that it be- 
haves like a single body; its equivalent Ji gravi- 
tational moment is negligible. Indeed, at distance 
r = 10 AU 



■h 



2 \r ) (mi + m^) 2 



10" 



(2) 



This can be confirmed easily by a direct numerical 
integration. The difference between trajectories 
computed for three-body (1, 2, 3) and two-body 
(1+2, 3) configurations is insignificant and always 
smaller than observational uncertainties (see Fig- 
ure 

We also make use of the following two con- 
straints: (i) initial positions 2^3, yhs and zero 
time To of the 3rd body correspond to a selected 
speckle-interferometry datum (e.g., the mean of 
the first two points, or to the third point); (ii) 3rd 
body initial velocity components are almost tan- 
gent to the observed interferometric trajectory in 
the (x, y) plane. 

Initial conditions of the integration are speci- 
fied in an arbitrary (usually 1+2-centric) frame. 
We then perform a transformation to a barycen- 
tric frame. The numerical integration runs in the 
barycentric Cartesian frame, where x, y axes cor- 
respond to sky plane, z axis is oriented from the 
observer towards the system. We use AU, AU/day 



units for positions and velocities. 

We integrate the system forward for 10,000 days, 
and backward (i.e., with opposite sign of initial 
velocities) for 20,000 days in order to cover the 
observational time span. The time step used is 
At = 10 days and the precision parameter of the 
BS integrator is e = 10~ 8 . Finally, we transform 
the output back to the (l+2)-centric frame and 
linearly interpolate the output data to the exact 
times of observations. 

In order to compare the observations to our 
model we constructed a % 2 metric as follows: 



X 



<^sky 



Xlite ' 



V 2 
A rv ) 



(3) 



where 



Xsky 



cv « 3 -xhM) 2 + (y'h 



hZ 



yhMf 



(4) 

We denote x' h3 , y' h3 (l+2)-centric coordinates of 
the 3rd body calculated from our model, which 
were linearly interpolated to the times i s ky[*] of 
observations a;/^^], 2/^3 [«]. Distance d is used to 
convert angular coordinates to AU. Secondly, 



Mite 

Xlitc = 

i=l 



( Z bl+2 - Z bl+2\i})' 



(5) 



'litcL 



where z' bl+2 are barycentric coordinates of the 
(1+2) body computed from our model and inter- 
polated to the times tuto [*] of observations 2514.2 [i]. 
Because of freedom in the period determination 
and freedom in the selection of initial velocities, 
we have to detrend the light-time effect data (by 
two LSM fits of z' bl+2 (t) and 25142 (t)). Finally, 



iV rv 

£ 



Kh3 ~ V zh 3 [i]Y 



(G) 



where we again interpolate our model to the times 
t IV [i\. Note that in case of a 4-body configuration 
we will attribute the velocities to the 4th body and 
change this metric correspondingly (see below). 

Optionally, we can add an artificial function 
to x 2 m order to constrain the mass 7714 within 
reasonable limits, e.g., 



7714 



m 4n 



TO 4n 



100 



777 4ll 



777 4n 



(7) 
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Table 3: Free parameters of our dynamical 4-body model. 



no. 


parameter 


brief description 


1. 


d 


distance of the V505 Sgr barycentre 


2. 




mass of the 3rd body 


o 
0. 


Zh3 


position, (l+2)-centric, epoch To 


4. 


V x h3 


velocities, (l+2)-centric, epoch T 


5. 


Vyh3 




6. 


Vzh3 




7. 


?Ti4 


mass of the 4th body 


8. 


XhA 


positions, (l+2)-centric, epoch T 


9. 


Vh4 




10. 


Zhi 




11. 


Vxhi 


velocities, (l+2)-centric, epoch T 


12. 


Vyhi 




13. 


Vzh4 





Table 4: Fixed (assumed) parameters of our model. 



no. 


parameter 


brief description 


14. 


TOl+2 = 3.4 Mq 


mass of the (1+2) body 


16. 


Xhi+2 = OAU 


positions of the (1+2) body, 


17. 


Vhl+2 = 


(l+2)-centric 


18. 


Zhl+2 = 




19. 


Vxhi+2 = AU/day 


velocities 


20. 


= 




21. 


Vzhl+2 = 




22. 


Xh3 


positions of the 3rd body, 


23. 


Uh3 


(l+2)-centric 


24. 


T = 2446282.24375 JD 


UTC time corresponding 




(or 2447607.5185) 


to initial conditions 




Fig. 1. — Comparison of two 3rd-body trajectories, computed for three-body (1, 2, 3) and two-body (1+2, 
3) configurations. Left: an overview of the trajectories in a 1-centric frame. Right: a detail of the small part 
of the trajectory, where the difference is visible. Error bars denote speckle-interferometry observations. 
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with m4 m ; n — 0.1 M Q , 

^'4 max = 1.2 M e . The up- 
per limit follows from the fact that no other bright 
star is observed in the vicinity of V505 Sgr. 

A similar expression can be used to constrain 
the absolute value of velocity V4 (e.g., to be smaller 
than the escape velocity from the system, other- 
wise, we often obtain hyperbolic velocities). 

Occasionally, we use a different metric instead 
of ©: 

X 2 = W sky X^ ky + WntoXnto + ^rvX™ ( 8 ) 

with weights u> s ky > wutc , Wrv , in order to fit 
the interferometric trajectory better. There are 
only five points after the periastron passage, which 
would otherwise have too low statistical signifi- 
cance compared to a lot of light-time data. 

What can we expect about the 13-dimcnsional 
function x 2 (d, m 3l z h3 , . . . , v zhi )7 It will surely 
have many local minima, which would be sta- 
tistically almost equivalent. (One can shoot the 
4th body from a slightly different position with a 
slightly different velocity to get almost the same 
result.) The problem is degenerate in this sense. 
Clearly, there are strong correlations, e.g., be- 
tween the mass 1114 and the minimal distance of 
a close encounter (and consequently initial posi- 
tions/velocities of the 4th body). Minimisation of 
the x 2 function is thus a difficult task. 

We use a simplex algorithm (Press et al. 1997) 
to save computational resources and to find lo- 
cal minima. However, it is not our goal to find a 
global minimum of x 2 ; because of the degeneracy 
and the immense size of the parameter space. We 
anyway do not expect a deep, statistically signif- 
icant global minimum. Instead, we will choose a 
set of starting points for the 4th body and look for 
a subset of allowed solutions. 

On the other hand, in case we test a 3-body 
configuration only, the problem is much simpler: 
the six-dimensional x 2 (d, m 3 , z h3 , v xh3 , v yh3 , v zh3 ) 
is well-behaved and we may expect to find a unique 
solution (and its uncertainty). 

4. Results 

In the following subsections we consider and 
analyse several hypotheses about the nature of the 
V505 Sgr system: (i) there are three bodies only 
in V505 Sgr; (ii) the 3rd body directly perturbs 
the central pair; (iii) a steady mass transfer causes 



minima timing variations; (iv) there is modulation 
of mass transfer by the 3rd body; (v) a sudden 
mass transfer occurred around 2000; (vi) Appel- 
gate's mechanism is operating; (vii) a 4th body is 
present (either on a bound or hyperbolic orbit). 

4.1. The 3rd body alone on a Keplerian 
orbit 

At first, let us test a standard "null" hypothesis, 
i.e., only 3rd body exists (1714 = 0). It is possible 
to fit speckle data alone (loute = w TV = 0) by an el- 
liptical orbit with a (29±l)-yr period, especially, if 
we assume the first two 1985 measurements are er- 
roneous (offset by 50 mas, see Figure [3J left). The 
X 2 ky = 50 for this fit and the respective number 
of data points is iV s k y = 20. (Thought ideally, x 2 
should be comparable to N.) 

Note the x 2 k y would be much higher, if we 
include the 1985 measurements: xLy = 210, 
iV s ky = 22. It means, if these two measurements 
are not systematic errors, the 29-yr Keplerian or- 
bit is essentially excluded! The two respective 
measurements were obtained by two different tele- 
scopes during two different nights (see McAlister 
et al. (1987a) and McAlister et al. (1987b)). We 
checked measurements of another 34 stars in these 
publications, observed with the same telescope 
and during the same night as V505 Sgr, and we 
have found no indication of a wrong plate scale — 
all measurements lie on Keplerian ellipses within 
usual observational uncertainties (5 mas). We thus 
belive the 1985 measurements are not erroneous 
and they should be included in the x 2 metric. 

Without additional (non-positional) data it is 
not possible to distinguish between different incli- 
nations — there are equivalent low-/ and high-/ 
solutions with almost the same x 2 — 50. Never- 
theless, every inclined orbit of the 3rd body has 
to cause a corresponding light-time effect, other- 
wise must be considered wrong! Even a slight 
/ > 2° inclination would be easily detectable in 
the light-time effect data (see Figure [2j middle). 
A period analysis of the O — C data (with Pe- 
riod04 program) also does not show a prominent 
29-yr period. On the other hand, there is a clear 
signal at P = 39 yr, with an amplitude of the peak 
A = 0.0092 d. 

If we assume the O — C data are indeed caused 
by a light-time effect, there is a strong disagree- 
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ment of the 29-yr Keplerian orbit with the light- 
time effect data (and also with radial velocities), 
even prior to 2000! If we try to fit the whole orbit 
and light-time effect data together, we would have 
Xsky = 1^7 and N s ^ y — 20, i.e., such an orbit is 
excluded with a high significance. There are also 
clear systematic departures between the observed 
interferometric data and calculated Keplerian or- 
bit. 

The only possibility is the inclination of the 3rd- 
body orbit is almost zero I < 2°, so we do not see 
any light-time effect at all. The observed O — C 
variations then must caused by an entirely differ- 
ent phenomenon (see next Sections 14.21 to 14.61 for 
a detailed discussion). 

Nevertheless, there still remains a strong dis- 
agreement with the observed high radial veloci- 
ties v ra d3 — fOkm/s, because a non-inclined orbit 
should have w ra d3 ^5 lkm/s. We have no solution 
for this problem (unless there is a 4th body present 
in the system, see Sections 14.71 to 14.91) . 

4.2. Direct perturbation of the 1+2 or- 
bital period by the 3rd body 

One may ask, if the observed variations in min- 
ima timings, which correspond to the changes of 
the period of the order |AP| ~ 10~ 5 d, could be 
caused by a direct gravitational perturbation of 
the tight central pair (1, 2) by the orbiting 3rd 
body. In periastron, the minimum distance is of 
the order ~ 10 AU. In order to test this possibil- 
ity, we use our dynamical model with three bodies 
1, 2 and 3 taken separately. A detection of minute 
changes of the orbital period requires a smaller 
time step and higher precision of the BS integrator 
(At = O.OId, e = 10" 12 ). The resulting osculat- 
ing orbital period changes during one periastron 
passage are shown in Figure [3] They are much 
smaller than AP < 10~ 7 d. An extremely close 
encounter (within less then 0.1 AU, which corre- 
sponds to 0.001 arcsec) would be needed to change 
the orbital period of the tight Algol system sub- 
stantially. 

Moreover, anything directly connected with the 
3rd body should conform to the 39 year period of 
the minima timings and this, according to Sec- 
tion 14.11 is in conflict with any 29-yr Keplerian 
orbit of the 3rd body. 



4.3. Effects of mass transfer between 1 and 
2 

Past photometric and spectroscopic observa- 
tions confirm the central pair of V505 Sgr is a 
classical semi-detached Algol system, with a less- 
massive secondary filling its Roche lobe (Chamb- 
liss et al. 1993). In case of a conservative mass 
transfer, the sum of masses is constant 

M t (t) + M 2 (t) = K , (9) 

as well as the orbital angular momentum 

A(t)M*(t)Ml(t) = C, (10) 

where A(t) denotes the actual separation of the 
stars. We can substitute current masses and sep- 
aration A = 7.1 Pq (Chambliss et al. 1993) into 
these equations, compute constants K, C and con- 
sequently the dependence A(Mi) (see also Fig- 
ure |4} 

A(Mi) = CM^ 2 (K - Mi)" 2 . (11) 

A smooth conservative mass transfer should in- 
crease orbital period steadily, since in the V505 
Sgr case the mass ratio has been reversed already 
(Mi > M2). On contrary, we observe an abrupt 
decrease of the period AP = —1.2 x 10~ 5 d after 
2000. We thus conclude a simple mass transfer 
cannot explain the observer minima timings. 

4.4. Modulation of mass transfer between 
1 and 2 during the 3rd body en- 
counter 

In this section, we test if the 3rd body is ca- 
pable to change the Roche potential of the cen- 

1992 1994 1996 1998 2000 2002 2004 2006 2008 



g 4e-07 - 




to 



2. -4e-07 - 

Q. | 1 1 1 1 1 1 1 

48000 49000 50000 51000 52000 53000 54000 55000 
JD - 2400000 

Fig. 3. — Simulated osculating orbital period P 
of the central binary (bodies 1 and 2), perturbed 
by the 3rd body. The periastron passage occurred 
in 2000 and the corresponding change of period 
is AP < 10~ 7 d. The observed values of |AP| ~ 
10~ 5 d are much larger than in this simulation. 
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Fig. 2. — Best-fit solution for the trajectory of the 3rd body, which corresponds to speckle interferometry 
data (but excluding 1985 measurements). Neither light-time effect nor radial velocities were fitted in this this 
case. Parameters of the 3rd body are: mj, = 1.17 M©, Xh3 = — 0.15 AU, ynz = — 25.9 AU, Zh3 = — 0.63 AU, 
v xh3 = 0.0037 AU/day, v yh3 = 0.0017 AU/day, v zh3 = 0.0000 AU/day for T = 2447607.5185 JD. The 
inclination of the orbit is very low in this case (/ = 1.5°). The resulting x^ky = 52, with the number of data 
points iVsky = 20. Note there is a strong disagreement of this Keplerian orbit with both O — C data and 
radial velocities (total \ 2 = 1700, N = 90). 



tral binary (bodies 1 and 2) in a such a way, that 
the mass transfer rate dM/dt (and consequently 
dP/di) changes by a substantial amount. We add 
a 3rd-body term to the Roche potential 

n(x,y,z) = ± + ^ + hl + q y 3 + q -^, (12) 
ri r 2 2 r 3rd 

where q = M2/M1 denotes the mass ratio and sim- 
ilarly (73rd = Ms/Mx. We see immediately, that 
relative change of the potential due to the 3rd 
body at distance r 3rd ~ 10 AU is 5Q/ft ~ 10~ 13 . 
We do not find it likely, that such a minuscule 
perturbation of the potential, and thus the re- 
lated tidal acceleration, could produce significant 
effects. Consequently, we cannot explain min- 
ima timings variations by the modulation of mass 
transfer. Finally, similarly as in Section 14.21 this 
effect would be also in conflict with a 29-yr Kep- 
lerian orbit of the 3rd body. 

4.5. A sudden mass transfer of Biermann 
& Hall (1973) 

According to Biermann & Hall (1973) a sudden 
mass transfer between the Algol components may 
result in a temporary decrease of the orbital pe- 
riod, even thought mass is flowing from the lighter 
component to the more massive. In our case, 
we would need dM/dt as high as ~ 10 _6 M©/yr 
to explain period changes |dP/df| ~ 10 _6 d/yr. 



Such a mass transfer rate seems to be too large 
compared to theoretical models (Harmanec 1970), 
dM/dt > 10~ 6 M©/yr are reached only during a 
very short interval of time, before the reversal of 
mass ratio. 

Another problem of this scenario is that we ob- 
serve rather smooth periodic variations of the min- 
ima timings before 2000, which do not seem to be 
entirely compatible with this mechanism, which 
may be more irregular in time. This phenomenon 
is also rarely confirmed by independent observa- 
tions. (It would require a very precise photometry 
on a long time scale, or a spectroscopic confirma- 
tion of circumstellar matter.) Today, this mecha- 
nism is not generally accepted as a major cause of 
minima timing variations among Algol-type sys- 
tems. 

4.6. Applegate (1992) magnetic mecha- 
nism 

Applegate (1992) proposed a gravitational 
quadrupole coupling of orbit and shape variations 
of a magnetically active subgiant (2nd compo- 
nent) can result in variations of the orbital period 
and hence minima timings. In this scenario, the 
observed 39-yr period would correspond to the 
period of the magnetic dynamo. 

The 2nd (G5 IV) star rotates quickly (1.2 d), 
it has a convective envelope in this evolutionary 



9 



stage and, presumably, there is a differential ro- 
tation and operating dynamo, which can result in 
a sufficiently strong magnetic field (10 4 G), nec- 
essary for Applegate's mechanism to work. Pe- 
riod changes of the order AP/P ~ 1CP 5 should 
also correspond to changes of the luminosity 
AL2/L2 = 0.1, in phase with minima timings. 
Unfortunately, we are not able to confirm this by 
our photometry (0.01 mag precision over tens of 
years would be required). 

In principle, this mechanism can explain min- 
ima timings variations, but it is not clear, why 
there is an abrupt change after 2000. An indepen- 
dent confirmation is rare and difficult. One of the 
possibilities might be a spectroscopic observation 
of magnetically active lines (Ca II H and K, or 
Mg II). This scenario also does not provide any 
solution for the observed large radial velocities. 

4.7. Distance, mass and the 3rd body or- 
bit (prior to 2000) 

Hereinafter, we assume minima timings varia- 
tions are caused mainly by the light-time effect 
due to the orbiting 3rd body. Because the orbit of 
the 3rd body prior the periastron passage in 2000 
seems unperturbed, we first determine the optimal 
distance d of the system, 3rd-body mass 7713 and 
orbit (2:3, v xh3l v y h3, v zh3 ). We use only the ob- 
servational data older than 2000 for this purpose. 

We compute x 2 values for the following set of 
initial conditions (we do not use a simplex here): 
d G (95,105)pc, Ad = lpc, m 3 G (1.1, 1.3) M Q , 
Am 3 = 0.1 M Q , z h3 G (2.0, 8.0) AU, Az h3 = 
1.0 AU, v xh3 £ (0.0033, 0.0040) AU/day, v yh3 G 
(0.0008, 0.0016) AU/day, v zh3 G (-0.0018,0.0012) 




0.5 1 1.5 2 2.5 3 3.4 



M, I M s 

Fig. 4. — The dependence of separation A of the 
eclipsing binary on the mass Mi of the 1st body, 
resulting from the conservation of total mass and 
orbital angular momentum. 



AU/day, Av xh3 = Av yh3 = Av zh3 = 0.0001 AU/day. 

The best-fit solution is displayed in Figure [5] 
Orbital period of the 3rd body is P = (39 ± 2) yr. 
The resulting distance d — (102 ±5) pc. This solu- 
tion is very similar to that in Mayer (1997). The 
parallactic distance of V505 Sgr given by Hippar- 
cos (tt = (8.40 ± 0.57) mas, d = 111 to 128 pc, cf., 
van Leeuwen 2007) is offset and even the error in- 
tervals do not overlap. 

Note that the radial velocities of the order 
— lOkm/s measured by Tomkin (1992) cannot be 
attributed to the 3rd body, which orbital veloc- 
ity should be much smaller ((—2.5 ± 0.5) km/s) 
according to interferometric and light-time effect 
data. Consequently, we do not fit the velocities 
in this case (w vv — 0), we are going to attribute 
them to the 4th body (in the next Section |4~8|) . 

Finally, it is important to mention that our so- 
lution does not depend on the two ("offset") 1985 
speckle measurements at all! We can exclude them 
completely from our considerations and the result 
would be the same. Our only assumption was that 
minima timings variations are caused by the light- 
time effect and this enforces the orbital period of 
P ~ 39 yr. (But coincidentally, both 1985 mea- 
surements fit perfectly this longer-period orbit.) 

4.8. Encounter with a 4th body (a \ 2 map) 

We next fix initial conditions of the 3rd body 
according to the results in Section 14.71 and model 
a perturbation by a 4th body under different ge- 
ometries. 

The free parameters of the model are: 771.4, 
xm, Vhi, ZhA, v xM , v y M, v zh4: - We include 
radial-velocity data, but we assume the spec- 
tral lines (and corresponding velocities) belong 
to the 4th body. We scan the following lim- 
ited set of initial conditions (over 8 million tri- 
als): 7714 G (0.5, 0.8) M Q , Am 4 = 0.05 Mq, x m G 
(38,45)AU, Ai M = 1.0 AU, y hi G (37,40)AU, 
Ay h4 = 0.5 AU, z M G (20,30)AU, Az hi = 
1.0 AU, v xM G (-0.011,-0.005) AU/day, v yM G 
(-0.010, -0.005) AU/day, v zhi G (-0.012, -0.006) 
AU/day, Av xhi = Av yhA = Av zhA = 0.0005 AU/day. 

A comparison of the best-fit solution with ob- 
servational data is displayed in Figure El We use a 
modified metric (jHJ with w s ^ y = 10, w\i tc = w TV = 
1. The respective trajectories of the bodies are 
shown in Figure [7] Note, however, that according 
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X/AU JD - 2400000 JD - 2400000 



Fig. 5. — Best-fit solution for the orbit of the 3rd body and the light-time effect before 2000: d — 102 pc, 
m 3 = 1.2 M , z h3 = 4.0 AU, v xh3 = 0.0036 AU/day, v yh3 = 0.0013 AU/day, v zh3 = -0.0015 AU/day for 
T = 2446282.24375 JD. The resulting x 2 = 88, with the total number of data points N = 46. The red lines 
denote differences between observed and calculated data. Radial velocities were not htted in this case, they 
are shown for comparison only. 



to the x 2 map (FigureHJ) there are many local min- 
ima, which cannot be distinguished from a statis- 
tical point of view, because the values of x 2 differ 
only little (x 2 G [284,325]). The corresponding x 2 
probabilities Q(x 2 |A), that the observed value of 
X 2 = 340 (for a given number of degrees of free- 
dom N — 105) is that large by chance even for 
a correct model, are too low (essentially zero). It 
may also indicate that real uncertainties might be 
a bit larger (by a factor of 2) than the values es- 
timated by us. Nevertheless, we will find better 
solutions using a simplex method (in Section |4~9| . 



4.9. 



Encounter with a 4th body (different 
geometry, simplex) 



We selected a different set of initial conditions 
for the following modelling. They serve as starting 
points for the simplex algorithm: 1714 = 0.5 Mq, 
x h4 G (-100,-10)AU, y M G (-50.1, -0.1) AU, 
z M G (0,50)AU, Ax M = Ay M = Az M = 
5.0 AU, v xM G (0.005, 0.015) AU/day, Av xhi = 
0.001 AU/day, v yhA G (0, 0.01) AU/day, Av yM = 
0.002 AU/day, v zhi G (-0.007, 0) AU/day, Av zhi = 
0.001 AU/day. The total number of trials reaches 
10 6 . 

We reject radial velocity constraints (w TV = 
0), although we can find a lot of allowed solu- 
tions with velocities in the correct range (v z h4 = 
—0.008 AU/day). On the other hand, we use a 
mass limit according to Eq. ([7])- An example of 
a typical good fit is shown in Figure |H1 We se- 




40^0 



Fig. 7. — Orbits of bodies (1+2), 3, 4 in a barycen- 
tric frame and their projection to the plane (z = 
-40 AU). 




Fig. 8. — Minima of the 7-dimensional function 
X 2 {rn i ,x hi ,Vhi,z hil v xM , v yhi ,v zh4 ) at given po- 
sitions (xhi,yh4:)- The minimum is taken over re- 
maining free parameters 7714, 27,4, i^m, Vj/M, v zh4- 
Cross is an overall minimum, black dots represent 
computed points. 
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X/AU JD-2400000 JD-2400000 

Fig. 6. — Best-fit solution for the trajectory of the 4th body, which best explains the observed trajectory 
of the 3rd body, light-time effect and radial velocities: 777,4 = 0.8 Mq, Xh4 — 45.0 AU, yh4 = 39.5 AU, 
z M = 28.0 AU, v xM = -0.0105 AU/day, v yM = -0.008 AU/day, v zM = -0.0075 AU/day. The resulting 
X 2 = 331, with the total number of data points N — 102. The motion of the 4th body captured in the 
left panel spans from 1994 to 2003. The squares connected by a straight line indicate the closest encounter 
between the 3rd and 4th body. 



lected one with mass around 7774 = 0.6 Mq; the 
corresponding \ 2 = 168, N = 73 and probability 
Q(X 2 \N) ~ 10" 9 , still too low. This solution can 
be further improved by a 15-dimensional simplex 
(i.e., with all parameters of the 3rd body free) to 
reach x 2 as l° w as 130 an d Q(x 2 \N) as high as 
10" 5 . 

As before, there are many solutions, which are 
statistically equivalent. We present allowed solu- 
tions in Figure[TU]as plots x 2 versus a free parame- 
ter, with each dot representing one local minimum 
found by simplex. Prominent concentrations of so- 
lutions in these plots can be regarded as an indica- 
tion of more probable solutions. Only minority of 
trials were successful. Most of them were stopped 
too early (at high x 2 ) due to numerous local min- 
ima. 

According to the histogram of masses 777,4 (Fig- 
ure [TTJ left) the values 777,4 < 0.5 Mq are less 
probable and the histogram peaks around 7774 = 
0.9 M©. Note the simplex sometimes tends to 
'drift' to zero or large masses, which leads to ar- 
tificial peaks at the limits of the allowed interval. 
The same applies to velocity V4. 

Histogram of total energies E4 of the 4th body 
(Figure [TTJ middle) shows a strong preference for 
hyperbolic orbits (E4 > 0), but elliptic orbits 
(E4 < 0) also exist (with a 1 % probability and 
slightly larger best x 2 = 199). The reason for this 
preference stems from the fact that 3rd body or- 



bit seems almost unperturbed prior to 2000, so one 
needs rather a higher- velocity encounter of the 4th 
body from larger initial distance. 

Typical minimum distances between the 4th 
and 3rd body during an encounter are around 
d m i n 3 ~ 6 AU and they are even smaller between 
the 4th body and (1+2) body d min i +2 ~ 1.5 AU 
(Figure [TT1 right). They are of comparable size 
and consequently a simple impulse approximation, 
i.e., an instantaneous change of orbital velocity, 
cannot be used to link the two elliptic orbits of 
the 3rd body (before and after the perturbation). 
There are no good solutions (with x 2 < 300), 
which would lead to an escape of the 3rd body. 

4.10. Observational limits of interferome- 
try and CFHT imaging 

Postulating an existence of a 4th body, inferred 
from its gravitational influence on the V505 Sgr 
system, we should check if this object could have 
been directly observed in the past. 

According to A. Tokovinin (personal commu- 
nication) the limit of recent interferometric mea- 
surements can be approximated by a linear de- 
pendence of brightness difference Sy in Stromgren 
y magnitudes on angular separation d of the com- 
ponents: Ay = 4.7mag at d = 0.15arcsec and 
Ay = 6.5 mag at d = 1 arcsec Adaptive optics at 
CFHT can reach even fainter. The limit in K 
band is given by Rucinski et al. (2007), Fig. 7 as 
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Fig. 9. — A typical solution (out of many) for the trajectory and light-time effect (without radial velocities): 
to 4 = 0.576 M Q , x m = 27.912 AU, y hi = -84.809 AU, z h4 = 30.482 AU, v xM = -0.00584 AU/day, v yhi = 
0.01612 AU/day, v z h4 — —0.00620 AU/day. The corresponding x 2 = 168, with the number of data points 
N = 73. This solution can be further improved by a 15-dimensional simplex to reach \ 2 as low as 130. 
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Fig. 11. — Left: Histogram of masses 1114 of all trials and good ones (with \ 2 < 300 and lower absolute 
value of velocity \vi\ < lOOkm/s). Middle: Histogram of total energies £4 of the 4th body at the epoch Tq. 
There is a strong preference for hyperbolic orbits (E4 > 0) , but one can find approximately 1 % of orbits 
with negative energies. Right: Histogram of minimum distances rf m ini+2, <inin3 of the 4th body from the 
(1+2) and 3rd body (only good trials are shown). According to these distributions, it is more likely, that 
the 4th body approached the (1+2) body closer than the 3rd body. Median values are <i m ini+2 = 1.6 AU 
and e? m i n 3 = 5.9 AU. 



a non-linear dependence AK(d). 

We can easily select solutions from Section |4~91 
which fulfil both limits, albeit a lot of them is ex- 
cluded by the CFHT limit (see Figure [12]). Note 
we are not able to predict exact magnitudes or 
positions of the 4th body, because there are still 
many solutions possible. 

Note there is an object in the USNO-A2.0 cat- 
alogue, very close to V505 Sgr: 0750-19281506, 
RAj 200 o = 298.277034°, DE J20 oo = -14.603839°. 
This corresponds to an angular separation 2.6 arcsec 
and position angle 235° with respect to V505 Sgr, 
at the epoch of observation 1951.574. The magni- 
tudes R = 10.9 mag and B = 11.8 mag are marked 
as uncertain (since the object is located in the area 
flooded by light of V505 Sgr). This is an inter- 
esting coincidence with "our" 4th body, but we 
doubt the source is real. Moreover, if the bright- 
ness of the USNO source is correct within ±1 mag, 
it should be above the observational limits. 

4.11. Constraints from spectral lines radial- 
velocity measurements 

In previous Section 14.81 we tried to attribute 
the observed high radial velocities to a hypothetic 
4th component. We thus have to ask a question: 
could the low-mass 4th component be visible in 
the spectrum? 



To this end, we used a grid of synthetic spec- 
tra based on Kurucz model atmospheres, which 
was calculated and provided for general use by 
Dr. J. Kubat (for details of the calculations, cf., 
e.g., Harmanec et al. 1997). We calculate syn- 
thetic spectra for 3 and 4 lights (stars) and com- 
pare them with the spectrum observed by Tomkin 
(1992), Fig. 2. This spectrum was taken at HJD 
= 2444862.588, close to the primary eclipse of the 
central binary, which decreases the luminosity of 
the 1st component and thus weak narrow lines of 
the 3rd (or 4th) component are more prominent. 

Modelling of spectra (relative intensities) re- 
quires a number of parameters: luminosities, ef- 
fective temperatures, surface gravity, rotational 
and radial velocities. Luminosities of the known 
components (out of eclipse) are: L\ = 26 Lq, 
L 2 = 3.8 Lq, L 3 = 2.1 Lq. The amplitude of 
the lightcurve is Am = 1.1 mag (Chambliss et 
al. 1993). The effective temperatures are ap- 
proximately (Popper 1980): T cS1 ~ 9000 K (cor- 
responding to A2 V spectral type) , T e g 2 — 6000 K 
(F8 IV to G6-8 IV), T cB3 ~ 6000 K (F8 V). We 
assume the following values of the surface gravi- 
tational acceleration: log gi — \ogg 2 — 4.0 (cgs 
units), log <?3 = 4.5 (valid for stars close to the 
main sequence). Rotational velocities of the 1st 
and 2nd components, a semi-contact binary with 
an orbital period 1.2 day, are synchronised by tidal 
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Fig. 12. — Observational limits (black lines) of interferometric measurements (left panel) and the CFHT 
imaging (right panel). Our solutions from Section [4791 are plotted as dots. Many of them are well below the 
observational limits. 



lock and are of the order v ro ti — v Iot 2 — lOOkm/s. 
These are in concert with the observed width of 
broad spectral lines AA = 6 A. For the 3rd compo- 
nent, we assume a lower velocity u r ot3 = 20km/s, 
usual for main-sequence stars. This matches the 
width of sharp lines. Radial velocities of the 1st 
and 2nd components are close to zero because of 
the eclipse proximity (i> r adi = Vmd2 = 0). 

We assume the following reasonable parameters 
for the 4th component: T 4 = 4000 K or 5700 K, 
log ^4 = 4.0, v Iot 4 = 20km/s. We construct a x 2 
metric 

X = g C bs [*] 2 ' (13) 

where iobsl*] denote observed relative intensities, 
fobsH associated uncertainties and I' is a sum of 
synthetic intensities weighted by luminosities 

, _ £j=i I'(T cSj , logffj , Vratj) ■ Lj 

and of course Doppler shifted due to radial veloci- 
ties (A' = A bs[«](l — Vradj/c)) and interpolated to 
the required wavelengths A' using Hermite poly- 
nomials (Hill 1982). We use a simple eclipse mod- 
elling: we decrease L\ according to the Pogson 
equation to get the observed total magnitude in- 
crease Am. Errors er b s [*] were estimated from 
the scatter in small continua, a — 0.01. Artifi- 
cially small errors a — 0.003 were assigned to the 
measurements in the cores of the narrow lines, in 
order to match precisely their depths. 



We constructed a simplex algorithm (Press et 
al. 1997) with the following free parameters: L4, 
^rad3; ^rad4- Other luminosities and radial veloc- 
ities remain fixed. This simplex is well-behaved 
and converges to final values almost regardless of 
starting point. There is no reasonable improve- 
ment, if we let all 8 parameters (Lj, i> ra dj) to be 
free. 

The results for two different temperatures are 
shown in Figure [TJ The best fit for T 4 = 4000 if 
is L4 = 0.22 Lq, and it is marginally better that 
the fit with 3 lights only (i.e., with fixed L4 = 0). 
The luminosity corresponds roughly to the mass 
TO4 / Mq oc (L4/L0) 1 / 4 = 0.68, which seems rea- 
sonable with respect to the results in Section |4"791 

Note we used w ro t3 = 20 km/s for rotational ve- 
locity of the 3rd body. No reasonable solution was 
found for w ro t3 as high as 100 km/s, which would 
cause a strong rotational broadening and almost a 
'disappearance' of spectral lines of the 3rd body. 
It means, that a low-mass 4th body alone cannot 
produce deep sharp lines. We thus suspect, there 
is a blend of lines in the spectrum of observed by 
Tomkin (1992), which may originate on the 3rd 
and 4th body, with low and high radial velocities. 
However, observations with high spectral resolu- 
tion would be needed to resolve such blending. 

4.12. Constraints from the stellar evolu- 
tion of the eclipsing binary 

To assess the long-term evolution of V505 Sgr, 
we need some information about the age of the sys- 
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Fig. 13. — A synthetic spectrum with intensity 
normalised to continuum for 4 and 3 lights (thick 
and thin lines) and for the temperature of the 
4th body equal to 4000 K, compared with the ob- 
served spectrum (gray line, which thickness corre- 
sponds to observational uncertainties), taken from 
Tomkin (1992). The wavelength range, which we 
fitted, was 5580 A to 5607 A. The best solution for 
the luminosity of the 4th body is L4 = 0.22 Lq 
(with x 2 = 465, N = 239) The fit with 3 lights 
only is worse, with \ 2 value equal to 542. For 
T 4 = 5700 if we would have L 4 = 0.43 Lq and 
X 2 = 495. 



tern. An upper limit for the age can be estimated 
easily from masses of stars. The semi-detached 
central binary (bodies 1 and 2) has a total mass 
(3.4 ± 0.1) M Q . In order to evolve into the cur- 
rent stage, when the 2nd lighter component fills 
its Roche lobe, the original mass of the 2nd star 
had to be at least slightly larger than half of the 
total mass, i.e., M2 > 1.7 Mq. The evolution of 
radius is shown in Figure 114) we are mainly con- 
cerned with the large increase of radius, when the 
star leaves main sequence. Given the uncertain- 
ties of the masses and unknown metallicities, the 
upper limit for the age is (2.0 ± 0.5) Gyr. 

In order to find a lower limit, we have to check 
a minimum separation of the components first 
(cf. Eq. [TT] and Figure [4j. A minimum separa- 
tion occurs when Mi = 0.5K , in our case A m i n = 
(5.9 ± 0.1) Rq. This value is larger than the ra- 
dius of a 3.4 Mq star during the whole evolution 
on the main sequence. Thus the mass transfer had 
to start later, in the red-giant phase. 

The maximum mass of the 2nd star had to be 
slightly below the total mass, i.e., Mi < 3.4 Mq. 
According to the R(t) dependence (Figure Hljl . 
the red-giant phase starts at the age of (0.26 ± 
0.03) Gyr, which could be considered as a lower 
limit for the age of the V505 Sgr system. 

5. Conclusions 

Generally speaking, we are able to explain the 
observed orbit of the 3rd body together minima 
timings and radial velocities by a low-mass 4th 
body, which encounters the observed three bod- 
ies with a suitable geometry. There is no unique 
solution, but rather a set of allowed solutions for 
the trajectory of the hypothetic 4th body. It is 
quite difficult to find a solution for both speckle- 
interfcromctry and light-time effect data. There 
are a few systematic discrepancies at the 2er level, 
which cause the likelihood of the hypothesis to 
be low. Possibly, realistic uncertainties er s k y are 
slightly larger (by a factor 1.5) than the errors es- 
timated by us. 

Of course, there are other hypotheses, which 
do not need a 4th body at all (a sudden mass 
transfer, Applegate's mechanism, etc.), but none 
of them provides a unified solution for all obser- 
vational data we have for V505 Sgr. 

Further observations of the light time-effect 
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Fig. 14.— Evolution of radii for stars with different masses (M = 1.6, 1.7, 1.8 and 3.3, 3.4, 3.5 M Q , Z = 0.03) 
and metallicities (M = 1.7 and 3.4 M , Z = 0.01,0.02). This computation was performed by the program 
EZ (Paxton 2004). The upper limit for the age of the 2nd component is then 1.5 to 2.5 Gyr, while the lower 
limit is between 0.24 and 0.29 Gyr. 



during the next decade can significantly constrain 
the model. A new determination of the systemic 
velocity of V505 Sgr may confirm, that the change 
in the O — C data after 2000 resulted from an ex- 
ternal perturbation. (Tomkin's (1992) value was 
(1.9 ± 1.4) km/s.) Spectroscopic measurements of 
the indicative sharp lines would be also very help- 
ful to resolve the problem with radial velocities 
mentioned in the text. 

If we indeed observe the V505 Sgr system by 
chance during the phase of a close encounter with 
a 4th star, we can imagine several scenarios for its 
origin: 

1. A random passing star approaching V505 
Sgr on a hyperbolic orbit. The problem of 
this scenario is a very low number density of 
stars. If we take the value n* ~ 0.073 pc -3 
from the solar vicinity (Fernandez 2005), the 
mean velocity with respect to other stars of 
the order v rc i ~ 10 km/s and the required 
minimum distance of the order d ~ 10 2 AU, 
we end up with a mean time between two 
encounters r ~ l/(n±v Ic id 2 ) ~ 10 12 yr, thus 
an extremely unlikely event. 

2. A loosely bound star on a highly eccentric 
orbit, with the same age as other three com- 
ponents of V505 Sgr. Unfortunately, there 



is a large number of revolutions and encoun- 
ters (10 2 to 10 5 ) over the estimated age of 
V505 Sgr and the system practically cannot 
remain stable over this time scale (Valtonen 
& Mikkola 1991). 

3. A more tightly bound star on a lower- 
eccentricity orbit, which experienced some 
sort of a late instability, induced by long- 
term evolution due to galactic tides, distant 
passing stars, which shifted an initially sta- 
ble configuration into an unstable state, e.g., 
driven by mutual gravitational resonances 
between components. The problem in this 
case is that tightly bound orbits of the 4th 
body are very rare in our simulations, thus 
seem improbable. 

None of the scenarios is satisfactory. Neverthe- 
less, we find the 4th-body hypothesis the only one 
which is able to explain all available observations. 
Clearly, more observations and theoretical effort is 
needed to better understand the V505 Sagittarii 
system. 
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